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Abstract 

Approximate equations are derived for the motion of a gyroscope on the earth's gravitational 
field using the Einstein, Infeld, Hoffmann surface integral method. This method does not require 
a knowledge the energy-momentum-stress tensor associated with the gyroscope and uses only its 
exterior field for its characterization. The resulting equations of motion differ from those of previous 
derivations. 



Already by 1927 Einstein had realized that the field equations of general relativity con- 
tained information about the motion of the sources of the gravitational field. 11 A complete 



working out of this possibility was achieved by Einstein, Infeld and Hoffmann^] (EJH ) and 
elaborated later in two papers by Einstein and Infeldj^]. Their method used only the source- 
free field equations of general relativity and avoided any consideration of the interiors of the 
sources, assuming only that they are compact, i.e. their sizes are small compared to the 
distance between them and that their velocities are small compared to the speed of light c. 
The nature of the sources is characterized by their exterior fields alone. The main drawback 
to the EIH procedure is that it requires a large amount of tedious calculation. 

An alternate approach to deriving equations of motion was introduced by Fock 4] and later 
developed further by Papapetrou [sj]. This second method takes into account the interiors of 
the sources and makes use of the conservation laws 

T^, v (1) 

where T^ v is the matter tensor associated with the source. (Here and in what follows Greek 
indices take the values 0, 1, 2, 3, Latin indices take the values 1,2,3, the Einstein summation 
convention is assumed, a comma denotes ordinary differentiation and a semicolon denotes 
covariant differentiation.) While this method requires somewhat less calculations, it has 

n 

the drawback that one must specify T^ v . Papapetrou [q] later used the Fock approach to 
derive equations of motion for 'test' bodies moving in an external gravitational field and 
applied it to the case of a spinning test body. However, the equations obtained by him were 
insufficient to determine completely the dynamical variables X^ 1 and S^, the position and 
spin angular momentum of the body, appearing in them. As a consequence it was necessary 
to impose restrictions from the outside on the components of S^ u in order to close the system 
of equations. These restrictions all require the vanishing of its space-time components (the 
iO components) but differ in which frame of reference these components vanish. Various 
authors 0] have proposed different restrictions leading to different equations of motion for 
S^ u . In a completely different approach using the Principle of Equivalence, Weinberg [§!] has 
obtained yet another set of equations for this variable that differs from those based on the 
Papapetrou equations. It should be pointed out that neither of these derivations make full 
use of the Einstein field equations - Papapetrou used them only to derive the conservation 
laws exhibited in equation (JTJ) and Weinberg did not use them at all. A resolution of these 
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different results has taken on a certain urgency with the launch of the NASA-Standard 
Gravity Probe B (GPB) designed to measure the change in orientation of a small gyroscope 
in earth orbit. This experiment was first suggested by Schiff[l| in 1960, whose calculations 
were based on the spin equations of motion obtained using the supplementary conditions 
proposed by Corinaldesi and Papapetrou (CP). 

In an attempt to resolve these issues I have undertaken to derive approximate equations 
of motion for compact spinning bodies using methods similar to the ones used by EIH. Today 
two developments have made it possible to perform the calculations needed for this purpose 
without undo labor - hig h speed personal computers and symbolic manipulation programs 
such as Mathematica 10] and Maple Q] together with the wonderful program grtensorQ]. 
With these tools it is now possible to obtain equations of motion in a matter of minutes 
with the assurance of correctness that would have required days or even weeks to perform 
by hand. 

Since the equations of motion obtained using the EIH approach are only approximate, 
it is necessary to specify the system one wishes to apply these equations to, in this case a 
compact gyroscope orbiting the earth, and to identify the small dimensionless parameters 
associated with this system that will be used in the expansions employed. The gyroscope 
used in GPB ( there are actually four of them) consists of an almost perfect sphere of 
fused quartz with a radius r g = .019 m and a mass m = .075 kg. and an initial angular 
velocity uj = 27000 rad/s. The gyroscope was launched into a near perfect circular polar 
orbit (eccentricity = .0014) of radius R = 7027 km. These values will determine the relative 
importance of the terms in the approximate expression for the gravitational field to be used 
to evaluate the surface integral terms that arise in the EIH approach. 

The form of the field equations for the gravitational field g^ v to be used here are due to 
Landau and Lifshitz[13]. Exterior to the field sources they have the form 

U^ p , p = -gt LL ^, (2) 

where 

jjw = = (l/l 6 7r){-^V CT - <? W <n},a , (3) 

g = det(g^ u ) and tj J jJ a/ is the Landau-Lifshitz pseudotensor. Because of the antisymmetry 
of U pup in its last two indices, it follows that [/ Mrs s is a three-dimensional curl and therefore 
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when equation (J2J) is integrated over a two-surface in a t = constant hypersurface, one gets 

§ s {U^° fl + gt LL ^)n r dS , (4) 
where n r is a unit surface normal. In a like manner one gets 

$ s {{xW r °) >0 - (ZW°),0 + ^Li" - J/s'W 

+ (l/167r){^(^V° - g» r g u0 ), }n r dS = 0. (5) 

It is these two last equations that are used in the EIH procedure to obtain equations of 
motion. 

In order to use equations (jJJ) and (j3J) it is necessary to obtain solutions of the field 
equations corresponding to the type of system being considered. Since in the case of GPB 
no such exact solution exists it is necessary to use approximate ones. In the case of GPB 
the system consists of two bodies, the earth and the gyroscope. Since the ratio of the 
masses M and m of the earth and the gyroscope is 1.25xl0 -26 it is clear that we can 
ignore completely any effect the gyroscope has on the earth's motion. We can therefore 
take the earth to be at rest at the origin of an inertial frame characterized by coordinates 
{ct, x, y, z}. At the location of the gyroscope the earth's field has the dimensionless value 
MG/Rc 2 = 6.5 x 1CT 10 where G is the Newtonian gravitational constant and M is the mass 
of the earth. The gyroscope's contribution to the gravitational field consists of three parts: 
its static part mG/r g c 2 = 2.9 x 10 -27 , an induction part mGV/r g c 3 = 7.3 x 10~ 32 , where 
V = 7.5 x 10 3 m/s is the orbital velocity of the gyroscope, and a spin contribution. The 
latter contribution depends on the gyroscope's spin angular momentum S = 0.29 kg m 2 /s 
and is given by SG/r 2 c 3 = 2.1 x 10~ 33 . From these numbers we can form the three small 
dimensionless parameters to be used in the construction of the approximate gravitational 
field. The first of these is the slowness parameter e = V/c = 2.5 x 10~ 5 . The second 
one is the ratio of the gyroscope's monopole field to that of the earth at the surface of 
the gyroscope, el = mR/Mr g = 4.6 x 10 -18 . Finally, the third small parameter is the 
ratio of the spin angular momentum of the gyroscope to its orbital angular momentum 
e2 = = S/mVR = 7.4 xlCT 11 . 

Taking these considerations into account, using units with G = c = 1 and measuring 
masses in units of M and lengths in units of R, the components of the gravitational field 
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g^ u have the form 



and 



M 77? 

goo = 1 - 2e 2 2e 2 el — - 2e 4 e2 li V l (6a) 

r rl 



A/f 77? 

g l3 = 8^-1 - 2e 2 - - 2e 2 el-) - e 4 e2( 7l V- + 7 ^) (6b) 



777- V' 

g i0 = 2e 3 el— - + e 3 e2 7i (6c) 
rl 



where 

and where r 2 = x 2 + t/ 2 + z 2 is the distance from the earth's center with coordinates {0, 0, 0} 
to the field point {x,y,z}, rl 2 = x\\ + xl\ + xl\ is the distance from the center of the 
gyroscope with coordinates{i?i, R2, R3} to this field point and is the antisymetric density 
with values +1 or -1 depending on whether ijk is an even or odd permutation of 123 and zero 
otherwise. Here only those terms that are needed to determine the lowest order equations of 
motion for the Sk have been included in these expressions for the approximate components 
of g^y. It is to be noted that the gyroscope monopole and dipole contributions to g^ are 
taken to have the same effective centers so that no supplementary conditions are needed to 
determine the time dependent of the Sk- 

To obtain the above expressions for the gravitational field one makes use of the field 
of a stationary spinning, spherically symmetric, body given by gio above. Since the GPB 
gyroscopes are moving in the earth's gravitational field, it is necessary to boost the static 
field to the velocity of the moving gyroscope. This boost is responsible for the terms in the 
expressions for g o and above that depend on Sj. In addition it is necessary that these 
fields satisfy the harmonic coordinate conditions to an accuracy that insures that the g^ 
are in fact approximate solutions of the Einstein field equations. In the present case this 
requirement is satisfied if 

W=9ST)r = O^el) (8) 
This will be the case provided that 

^,0 = 0(e 3 ) (9) 

Finallyit is necessary to discuss the dependence on time of the dynamical variables Ri and 
In their original paper, EIH introduced what they called the slow- motion approximation 



by assuming that the source coordinates depended on the time t through the combination 
et. Their procedure is equivalent to what is known today as a multiple time formalism 
and will be used in what follows. This being the case, condition (JUJ) can be satisfied if one 
assumes that 

Si = s 0i + es u (et) (10) 

where s^i is independent of et . (In higher orders of approximation soi will in general depend 
on e 3 t .) 

All that remains is to substitute the above expressions for into the surface integrals 
in equations (JH) and (0) and evaluate the integrals over a sphere surrounding the gyroscope. 
Most of the terms so obtained will depend on the radius of the sphere chosen and so must 
cancel as a consequence of the field equations ((21). Those terms that are independent of the 
sphere radius must vanish as a consequence of the motion of the sources, here the gyroscope, 
and hence are the desired equations of motion. Evaluating the surface integrals in equation 
E] to 0(e 4 el) yields the Newtonian equations of motion for a particle moving in the earth's 
gravitational field 

where a prime denotes differentiation with respect to et. 

To get equations for the spin variables Si it is necessary to evaluate the surface integrals in 
equation© to 0(e 6 e2), a task that would have been beyond my abilities to perform without 
the help of Mathematica and grtensor. One finds that these equations can be written in 
vector form as 

d Sl /d{et) = ^^{-(s -R')R+19(s -R)R+ 16s (R-R)} (12) 

These equation are to be compared to the ones obtained by Cornaldesi and Papapetrou 
given by 

ds/dt = 2— {--(s-R)R + (s • R)R + 2s(RR) - g^^O 3 ' R ) R > ( 13 ) 

and Weinberg's equation 

da/d(t) = ^{-2(s-R)R + (s • R)R - 2s(R-R)} . (14) 
R 

where a dot over a quantity denotes differentiation with respect to t. These latter two sets 
of equations are considered to be exact by their authors but can be solved approximately 



by substituting for s the expression given for it in equation (fTUj) . It is clear that they will 
give different results than those obtained from equations (fT2|) . In all of these equations the 
dependence of R on t is gotten by solving the Newtonian equations of motion 
In the case of a circular motion in the xy-plane we can take, in equation 

Ri = Rcos(eujt) and R2 = Rsin(eu>t). (15) 

where uo is the angular velocity of the gyroscope in its orbit. To find the secular change in 
s with time we can average Equation. (|T2|) over an orbital period T = 2-k/uj. If one takes So 
to lie in the plane of the orbit and the xy-axes are chosen so that soi = sq and S02 = 0, the 
resultant change Asi is given by 

Asn=0 (16) 

and 

Asi2 = 2ns M/R (17) 

so that the angular change A8 in the direction of s is given by 

Ad = 2nM/R. (18) 

A similar analysis yields a value AO = M/2R for the CP equations and 3M/2R for the 
Weinberg equations. Why the difference in the three results? In the case of the Papapetrou- 
Corinaldesi equations the authors made assumptions concerning the matter tensor T^ v that 
are not justified and Weinberg relied on the Principle of Equivalence and identified the space 
part of a four-vector with an axial three- vector. 

It is also possible to take account of the spin-spin interaction between the gyroscope and 
the earth's rotation about its axis. To do so it is first necessary to introduce a fourth small 
parameter e3 = Q/u = 6.8 x 10~ 2 , where Q is the angular velocity of the earth, into our 
expansions. The gravitational field of the earth's rotation can be taken account of by adding 
to the expression (jHcjl for g i0 a term 

g i0 = e 3 e3e ijk ^j^ (19) 

where Sk = IQk is the earth's angular momentum and I is its moment of inertia about 
its axis of rotation. With this addition to the gravitational field, the surface integral (0) 
introduces an additional term ts in the equation of motion (|T2|) for si given by 



Ts= 2% ^ 2(So X S) " 3(S ' R)(S ° X R} ' (20) 



For a circumpolar orbit with xy-axes now chosen so that S = {S, 0, 0}, assuming that s = 
{so cos(</?), So sm(ip), 0} and after averaging over an orbital period one finds an additional 
change in Sx given by 

with a corresponding angular change in the direction of s given by 

A0 = ^ S[n{ ^ (22) 

It is amusing to think that if this additional change in direction could be measured with 
enough accuracy one could use the result to determine / by assuming that general relativity 
was the correct theory of gravity. 
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